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摘 要 本 文 针对 微 通 道内 气 粒 间 流 动 传 热 过 程 开 展 数值 研究 , 所 建 模型 中 气体 处 理 为 可 压缩 / 变 物性 流体 , 并 在 气 固 交界 
处 采用 速度 滑 移 和 温度 跳跃 边界 条 件 以 考 虚 其 微 尺度 效应 . 在 数值 研究 基础 上 , 分 析 了 微 通道 受 限 空间 、 克 努 森 数 、 气体 流 


速 以 及 颗粒 表面 温度 对 微 通道 内 气 粒 换 热 的 影响 . 研究 结果 表明 , 受 限 空间 结构 将 强化 气 粒 间 换 热 过 程 , 颗粒 表面 平均 传 
热 努 赛 尔 数 随 微 通道 气体 流量 增 大 而 增 大 , 克 努 森 数 增 大 以 及 颖 粒 表 面 温度 升 高 都 将 导致 颗粒 表面 平均 传 热 努 赛 尔 数 减 小 . 
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Numerical Simulation of Gas-particle Flow and Heat Transfer in Slip flow Regime 
in a Microchannel 


ZHOU Jun LIU Zhen-Yu WU Hui-Ying 
(School of Mechanical Engineering, Shanghai Jiao Tong University, Shanghai 200240, China) 


Abstract In this paper, a numerical model considering the compressibility and variable properties 
of gas fluid was established to simulate the gas-particle flow and heat transfer in the microchannel. 
The velocity slip and temperature jump boundary conditions were implemented at the gas-particle 
interface in the numerical modeling. The influences of microchannel confined space, Knudsen number, 
gas velocity and particle surface temperature on the heat transfer between gas and particle were 
studied based on the numerical predictions. The results show that the heat transfer is enhanced due to 
the confined space of microchannel. Moreover, the average Nusselt number increases with increasing 
gas flow rate, while it decreases with the increases of Knudsen number and particle surface temperature. 
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近年 来 , 随 着 微 加 工 制造 技术 水 平日 益 提 高 , 各 
种 微型 化 系统 不 断 涌现 . 其 中 , 微 纳 卫星 已 成 为 未 来 
航天 器 发 展 的 新 方向 04, 微型 固体 能 源 系统 是 实现 
微 纳 卫 星 在 轨 姿 控 的 新 型 推进 系统 , 具有 集成 度 高 、 
运行 可 靠 等 特点 ,受到 研究 者 高 度 关注 Pl. 对 于 微 
型 固体 能 源 推进 系统 ， 全 面 掌握 微 广 度 喷 管内 气体 
-固体 颗粒 流动 换 热 特性 是 其 设计 开发 成 败 的 关键 
问题 

微 尺 度 气 体 - 颗粒 流动 传 热 研 究 尚 处 于 初步 
UrBr. Tsuji 四 研究 了 气 固 两 相 之 间 动 量 、 能 量 厅 
合 机 制 以 及 不 同 尺 度 气 固 两 相 流 的 处 理 方法 . Barber 
等 回采 用 数值 模拟 和 理论 分 析 方 法 研究 低 Reynolds 
(Re) 数 微 通 道内 固定 颗粒 表面 搜 力 系数 (Ca) 问题 ， 
研究 发 现 阻塞 率 (D/H) 对 颗粒 表面 Ca 具有 较 大 影 
响 , 同时 得 到 颗粒 表面 Ca 的 表达 式 , 但 研究 中 并 未 
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考虑 温度 跳跃 的 影响 ，Mohajer 等 [9 分 析 气 体 变 
物性 对 无 限 大 空间 内 三 维 球形 颗粒 与 周围 空气 传 热 
特性 的 影响 ， 研 究 表明 当 颗 粒 表面 与 气体 来 流 温 度 
之 差 小 于 200 K 时 , 颗粒 表面 速度 滑 移 效应 对 传 热 
的 影响 大 于 温度 跳跃 效应 带 来 的 影响 ,颗粒 表面 平 
均 传 热 努 赛 尔 数 (Nu) 随 克 努 森 数 (Km) 增 大 而 增 
大 ， 而 当 两 者 温度 差 大 于 600 K 时 具有 相反 的 结 
iE, 但 该 研究 并 未 涉及 气体 可 压缩 效应 。Kishore 和 
Ramteke M 研究 了 速度 滑 移 边界 条 件 下 ， 球 型 颗粒 
与 气体 之 间 的 对 流 换 热 问题 ， 但 该 研究 中 并 未 考虑 
温度 跳跃 效应 对 传 热 的 影响 。 

综 上 所 述 ， 目 前 针对 微 通 道 受 限 空 间 内 气 粒 流 
动 传 热 的 研究 还 较为 匮乏 ， 现 有 模型 并 未 充分 考虑 
微 尺度 条 件 下 气体 的 可 压缩 性 对 流动 传 热 影 响 ， 对 
气体 变 物 性 影响 的 研究 也 不 够 深入 。 本 文 将 数值 模 
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拟 微 通道 受 限 空间 内 气 粒 间 流 动 传 热 过 程 ， 分 析 气 
固 边 界 处 速度 与 温度 非 连 续 性 对 气体 -颗粒 间 换 热 
过 程 作用 规律 ,研究 微 尺度 气体 可 压缩 性 及 变 物性 
对 气 粒 流动 传 热 特性 的 影响 程度 ， 
1 数学 模型 

随 着 微型 器 件 特征 尺寸 减 小 ， 其 内 部 气体 可 压 
缩 性 及 稀薄 作用 将 不 可 忽略 ， 气 体 流 动 的 连续 性 假 
设 将 不 再 成 立 591, 通常 引入 Kn 数 判断 微 尺 度 气 体 
是 否 适用 连续 性 假设 条 件 001; 
入 
L 


其 中 , L 为 流动 特征 尺寸 (本 文 取 颗 粒 直径 Dp), A 
为 气体 分 子平 均 自由 程 . 根据 Kn 数 大 小 可 将 气体 
流动 划分 为 四 个 区 : 连续 介质 区 (Kn < 10-?)、 温 
度 跳跃 与 速度 滑 移 区 (10 < Kn < 1071), HEK 
(107! < Kn «10*) 以 及 自由 分 子 区 (Km > 101) W, 
本 文 研究 的 气 粒 两 相 流动 传 热 问 题 属于 温度 跳跃 与 
速度 滑 移 区 , 气体 主流 区 域 Navier-Stokes 方程 和 傅 
里 叶 定 律 依然 适用 , 但 在 气 固 交 界 处 速度 与 温度 将 
出 现 非 连 续 特 征 . 

针对 稳 态 二 维 可 压缩 无 内 热源 牛顿 流体 (不 考 


Kn = (1) 


虑 气体 体积 力 影 响 ), 其 主流 区 控制 方程 组 : 
连续 性 方程 : 
O(pu) (pv) — 
Ox T Oy | i (2) 
2 方向 动量 方程 : 
9 (pu?) O(puv) ^ Op Trs , OTy« 
-ar 1 Oy T Os! 0s ^ By O 
y 方向 动量 方程 : 
ð (pv?) ð (puv) E | Op 十 OTyy 十 Tay (4) 
Oy Ox Oy Oy ox 
能 量 方程 : 
D[IoEvu 十 pu | 8|pEv + pv] 
De SNC A 十 ——— — 
Oz Oy 
oT oT (5) 
o M F o [Euro ru 
Oz Oy 
| 2u /0u Ou .2u(j0v Ov 
RP, Tes = 3 e x) Tw 3 e x) i 


Te = tey = n(A E) u, 分别 表 示 气 体 在 


z, 8 方向 上 的 速度 分 量 , k 为 气体 导热 系数 , ENT 
) 


体 动能 (V?/2) 与 内 能 (e) 之 和 : 
V? 


微 通 道 进口 设 定 为 速度 入 口 边 界 条 件 ， 出口 采 
用 压力 出 口 边界 条 件 。 当 Kn < 0.1 时 ， 可 采用 
如 下 边界 条 件 描述 其 气 固 相 间 速 度 与 温度 不 连续 效 
JY. Maxwell 一 阶 速度 滑 移 边 界 条 件 : 
2 一 au 入 


Ug 一 Uw = id =T wan 
~ (7) 
3 —1 
FL Q ) (—qt) "m 
4 qypRTilip 


其 中 , ws 为 近 壁 面 气体 速度 , uw 为 壁面 速度 .cv 为 
切 向 动量 协调 系数 (在 本 文中 取 1), / 为 空气 动力 黏 
BE, 和 为 气体 分 子平 均 自由 程 . 式 (7) 等 式 右 侧 第 一 
MAAR EN, 第 二 项 为 热 蠕动 项 RERE 
项 相 比 热 蠕 动 项 可 以 忽略 97), 因此 , 本文 研 究 中 采 


用 的 速度 滑 移 边界 条 件 为 : 
= on (8) 
Ug — Uw = " n wall 


Smoluchowski 一 阶 温度 跳跃 边界 条 件 : 
2—or 2y A 0T (9) 

or y+1Pr ón 
其 中 , Te 为 近 壁 面 气体 的 温度 ，7Tw 为 壁面 温度 . 
a, ?分别 为 热 协 调 系数 和 气体 比热容 比 (本 文 研 究 
中 分 别 取 1 1.4). 

通过 模拟 平板 间 充 分 发 展 的 微 尺 度 Poiseuille 流 
动 来 验证 本 文 所 建 模 型 的 可 靠 性 . 如 图 1 所 示 , 将 
数值 模拟 结果 与 理论 结果 UU 进行 比较 ,数值 计算 
得 到 的 平板 间 无 量 纲 温度 和 无 量 纲 速度 分 布 与 理论 
解 结果 吻合 良好 . 
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(a) Dimensionless temperature distribution 
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(a) Dimensionless velocity distribution 


图 1 微 尺度 Poiseuille 流动 数值 解 与 理论 解 对 比 


Fig. 1 Comparation between numerical and analytical results 


of micro-scale Poiseuille flow 


2 物理 模型 


针对 本 文 所 研究 的 微 通 道 气 粒 流动 传 热 问 题 ， 
计算 物理 模型 如 图 2 所 示 . 于 30 umx10 um 的 微 通 
道中 心 位 置 男 定 一 个 高 温 圆 形 固体 颗粒 ， 常 温 气 体 
以 恒定 速度 从 左 侧 入 口 流 入 通道 冷却 该 颗粒 。 对 气 
体 流动 区 域 进行 网 格 划 分 ， 在 微 通 道 壁面 和 颗粒 表 
面 附近 区 域 采 用 贴 体 细密 网 格 以 提高 计算 精度 . HB 
先 对 不 同 Kn 数 流动 条 件 的 计算 网 格 进行 无 关 性 检 
Jy, 表 1 给 出 了 Kn—0.03 的 检验 结果 , 在 综合 考虑 
计算 速度 与 精度 基础 上 , 选择 网 格 3 进行 Kn—0.03 
条 件 下 的 流动 传 热 计算 ， 其 他 Kn 数 条 件 下 也 采用 


10 um 


L-30 um 
图 2 计算 物理 模型 示意 图 


Fig. 2 Schematic of computational model 


表 1 网 格 无 关 性 验证 (Kn—0.03) 
Table 1 Verification of grid-independency( Kn 


—0.03) 
网 格 总 数 ” ”平均 Nt ë MARA 
网 格 1 8875 1.1721 0.0949 
网 格 2 13212 1.0979 0.0257 
网 格 3 25034 1.0724 0.0019 
网 格 4 33655 1.0704 0.000 
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相同 的 步 又 进行 网 格 的 划分 。 在 此 基础 上 ， 数 值 模 
拟 了 Kn=0.01~0.1, 颗粒 表面 温度 Tp—500-1100 K 
以 及 微 通道 空气 进口 速度 Vi,=1~7 m/s 工 况 下 微 通 
道内 气 粒 间 的 流动 与 传 热 过 程 . 


3 结果 与 讨论 


通过 比较 微 通 道内 单个 颗粒 和 无 限 大 空间 内 单 
个 颗粒 与 气体 之 间 的 换 热 特性 ， 研究 微 通道 受 限 空 
间 对 颗粒 与 气体 之 间 换 热 特性 的 影响 . 图 3 给 出 了 
颗粒 V, —3 m/s, T,—900 K 工 况 下 , 颗粒 表面 平均 传 
热 Nu 数 随 Kn 数 的 变化 情况 . 从 图 3 可 以 看 出 , 在 
考虑 气体 变 物 性 条 件 下 ， 置 于 无 限 大 空间 中 单个 颗 
粒 表面 平均 传 热 Nu 数 随 Kn 的 增 大 而 减 小 . 同样 ， 
微 通道 受 限 空间 内 颗粒 表面 平均 传 热 Nu 数 也 随 Kn 
数 的 增 大 而 减 小 , 但 是 减 小 的 幅度 大 于 无 限 大 空间 
单个 颗粒 情况 .这 是 因为 在 Kn 数 较 小 (Dp 较 大 ) 
时 ， 受 限 空间 内 空气 流通 面积 较 小 ， 导 臻 颗粒 表面 
附近 空气 流速 较 大 , 因此 会 增强 流动 换 热 过 程 ; 而 随 
Kn BOE, 受 限 空 间 对 空气 流速 影响 减弱 , 逐渐 接 
近 无 限 大 空间 内 单个 颗粒 的 情况 ,因此 两 种 情况 的 
颗粒 表面 平均 传 热 Nu 数 也 逐渐 接近 . 图 4 给 出 了 
Kn—0.05, V, —3 m/s, Tp=900 K LU F, 无 限 大 空 
间 单 个 颗粒 周围 气体 速度 与 温度 分 布 . 图 5 给 出 了 
Kn—0.05, V, —3 m/s, Tp=900 K 工 况 下 , 微 通道 受 
限 空间 内 单个 颗粒 周围 气体 速度 与 温度 分 布 。 由 图 
4(a) 与 图 5(a) 可 以 看 出 ， 由 于 微 通道 壁面 的 存在 ， 
颗粒 周围 气体 在 受 限 空间 中 的 流速 明显 大 于 无 限 大 
空间 内 颗粒 周围 气体 的 流速 ， 而 空气 流速 的 增加 有 
利于 颗粒 与 其 周围 气体 间 的 换 热 过 程 . 从 图 4(b) 与 


通道 受 限 空间 
民 大 空间 


0.01 0.03 0.05 0.07 0.09 0.11 
Re 
图 3 受 限 空间 对 颗粒 表面 传 热 Nu 数 的 影响 (Vis —3 
m/s, Ty —900 K) 
Fig. 3 The influence of confined space on average Nu 
number(Vin=3 m/s, 7T, —900 K) 
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图 5(b) 可 观察 到 ， 由 于 受 限 空间 中 气体 流速 较 快 ， 
能 及 时 将 热量 带 走 ， 因 此 受 限 空 间 中 颗粒 周围 气体 
温度 低 于 无 限 大 空间 中 颗粒 周围 气体 温度 . 


Velocity: 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 55 6 
(a) 速度 分 布 
(a) Velocity distribution 


Temperature: 350 400 450 500 550 600 650 700 750 800 850 
(b) 温度 分 布 
(b) Temperature distribution 
图 4 无 限 大 空间 内 颗粒 周围 速度 与 温度 分 布 
(Kn=0.05, Vi; 3 m/s, Ty —900 K) 
Fig. 4 Velocity and temperature distribution around micro 
particle in infinite space (Kn—0.05, Vi, —3 m/s, Ty, —900 K) 


为 进一步 研究 颗粒 表面 平均 传 热 Nu 数 与 颗粒 
直径 (Dp) 间 的 关系 . 图 6 给 出 了 Tp=700 K 工 况 
下 , 颗粒 表面 平均 传 热 Nu 数 随 Kn 数 变化 情况 . 一 
方面 在 空气 变 物性 条 件 下 ， 温 度 跳跃 对 传 热 削弱 程 
度 大 于 速度 滑 移 增强 程度 ， 因 此 随 着 Kn 数 的 增加 
颗粒 表面 平均 传 热 Nu 数 减 小 。 另 一 方面 ， 由 于 本 
文中 计算 Kn 数 的 特征 尺寸 取 Dp, 因此 随 着 Kn 数 
不 断 增 大 (De 不 断 减 小 ), 微 通道 内 流体 流通 面积 
大 ， 相 同 微 通道 Vu 进口 条 件 下 ,颗粒 表面 气体 速 
度 随 之 减 小 ,因此 颗粒 表面 传 热 减弱 .综合 以 上 两 
种 因素 ，Kn 数 增加 将 引起 颗粒 表面 平均 传 热 Nu 数 
减 小 。 

微 通道 Vi, 会 直接 影响 流 经 颗粒 表面 气体 速度 ， 
而 颗粒 表面 气体 流速 与 其 传 热量 密切 相关 . 图 7 给 
出 了 Kn=0.03 和 Kn=0.05 两 种 情况 下 微 通 道内 颗 
粒 表 面 平 均 传 热 Nu 数 随 微 通 道 Vu 的 变化 趋势 . 从 


图 中 可 以 看 出 ， 在 相同 颗粒 表面 温度 条 件 下 ， 由 于 
进口 速度 的 增加 ， 流 经 颗粒 表面 气体 的 速度 也 随 之 
增加 ， 有 利于 换 热 . 因此 ， 随 微 通道 Vi, 不 断 增 加 ， 
颗粒 表面 平均 传 热 Nu 数 也 随 之 增加 。 在 图 7 中 还 
可 以 观察 到 ,7P 较 高 时 , 颗粒 表面 平均 传 热 Nu 数 
BE Vi, 增加 的 幅度 小 于 Te 较 低 的 情况 . 这 是 因为 Tp 
较 高 时 , 温度 跳跃 效应 对 传 热 的 削弱 作用 也 较 强 , 因 
此 增加 Vin 并 不 能 大 幅度 增强 换 热 ， 


Velocity: 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 


(a) 速度 分 布 
(a) Velocity distribution 


Temperature: 350 400 450 500 550 600 650 700 750 800 
(b) 温度 分 布 
(b) Temperature distribution 
图 5 微 通道 受 限 空间 内 颗粒 周围 速度 与 温度 分 布 
(Kn=0.05, Vi, 23 m/s, Ty —900 K) 


Fig. 5 Velocity and temperature distribution around micro 


particle in a confined microchannel (Kn=0.05, Vi, —3 
m/s, Ty —900 K) 
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图 6 Kn 数 对 颗粒 表面 传 热 Nu 数 的 影响 (Tp —700 K) 
Fig. 6 Influence of Kn number on average Nu number 
(Tp—700 K) 
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(a) Kn=0.03 
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Nu 


Vi (m/s) 
(b) Kn-0.05 
图 7 微 通 道 进口 速度 (Vin) 对 颗粒 表面 传 热 Nu 数 的 影响 


Fig. 7 Influence of inlet velocity on average Nu number 


气体 的 物性 对 颗粒 与 气体 之 间 传 热 有 重要 影响 ， 
因此 研究 与 空气 物性 密切 相关 的 颗粒 表面 温度 (Tp) 
对 颗粒 表面 传 热 的 影响 具有 重要 意义 . 图 8 给 出 了 
Kn-—0.03 和 Kn—0.05 两 种 情况 下 , 颗粒 表面 平均 传 
热 Nu 数 随 Tp 的 变化 情况 . 随 Tp 增加 , 颗粒 与 流 
经 颗粒 表面 空气 间 的 温差 增 大 ， 相 应 的 温度 梯度 随 
之 增 大 , 温度 跳 牙 效应 也 逐渐 增强 ， 传 热 热 阻 增 大 ， 
不 利于 换 热 ; 而 随 着 Te 增加 , 空气 的 导热 系数 也 相 
应 增 大 并 有 利于 换 热 . 上 述 两 者 存在 的 “竞争 ”的 关 
系 决定 了 最 终 的 传 热 效 果 。 由 于 温度 跳跃 效应 对 传 
热 的 削弱 作用 大 于 导热 系数 产生 的 增强 作用 , 因此， 
B& Tp 的 升 高 , 颗粒 表面 平均 换 热 Nu 数 减 小 . 同时 ， 
在 图 8 中 还 可 以 观察 到 , 在 Tp 较 低 时 , 由 于 温度 跳 
路 对 传 热 的 削弱 作用 有 限 ， 因 此 Vin 较 大 时 的 换 热 
强度 明显 大 于 Vin 较 小 的 换 热 强 度 , 而 随 Tp 逐渐 升 
高 时 ， 由 于 温度 跳跃 对 传 热 的 前 弱 作 用 增强 ， 各 个 
Vin 工 况 下 的 颗粒 表面 平均 传 热 Nu 数 逐 渐 接 近 ， 


—IK—V,-1 m/s —9—y,.—2 m/s 
一 全 一 V =3 m/s——V,,-4 m/s 
一 多 一 V-5 m/s—$— V;,—6 m/s 
—9— V,,—7 m/s 


500 600 700 800 900 1000 1100 


(a) Kn=0.03 
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图 8 颗粒 表面 温度 (TP) 对 颗粒 表面 传 热 Nu 数 的 影响 
Fig. 8 Influence of particle surface temperature on average 
Nu number 


4 t5 论 


本 文 针对 微 通道 内 气 粒 流动 传 热 过 程 开展 数值 
模拟 研究 , 所 建 模型 中 气体 处 理 为 可 压缩 / 变 物 性 流 
体 ， 并 在 颗粒 表面 以 及 微 通道 壁面 处 采用 速度 滑 移 
和 温度 跳 路 边界 条 件 以 考虑 其 滑 移 区 微 尺 度 效 应 
研究 表明 ， 受 限 空 间 对 颗粒 与 空气 间 传 热 有 明显 影 
响 ， 随 着 Kn 数 的 增 大 (De 的 减 小 ), 颗粒 表面 平均 
传 热 Nu 数 随 之 减 小 . 同时 , 颗粒 表面 平均 传 热 Nu 
数 随 微 通道 Vu, 增加 而 增 大 ， 随 颗粒 Te 增 大 而 减 
小 . 但 相 较 Dp, 微 通道 Vi, 及 Tp 对 颗粒 表面 平均 
传 热 Nu 数 影响 较 小 ， 
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